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We consider, within the statistical hadronization model, the near central rapidity y ~ integrated 
hadron yields expected at LHC -^snn = 2.76 TeV ion reactions, for which the total charged hadron 
rapidity most central head-on collision yield is dh/dy\y=o ~ 1800 . For the chemical equilibrium 
SHM, we discuss composition oidh/dy as function of hadronization temperature. For chemical non- 
equilibrium SHM, we input the specific strangeness yield s/ S, demand explosive disintegration and 
study the break up as a function of the critical hadronization pressure P. We develop observables 
distinguishing the hadronization models and conditions. 

PACS numbers: 12.38.Mh, 24.10.Pa, 25.75.-q 



Introduction: According to the theory of strong in- 
teractions, Quantum Chromo-Dynamics (QCD), quarks 
and gluons are confined inside hadrons. At sufficiently 
high temperature lattice computations demonstrate that 
the deconfined Quark-Gluon Plasma (QGP) prevails fH]. 
We seek to prove this QCD paradigm in high energy 
heavy ion collisions where heavy nuclei are crushed on 
each other, forming a small drop of thermalized decon- 
fined QGP matter. This drop hadronizes into a high 
multiplicity of particles and we seek to determine the 
physical properties of QGP considering this multiparti- 
cle final state. 

At LHC energy ^snn = 2.76 TeV, for near head-on 
5% most central Pb-Pb collisions, the pseudo-rapidity 
density of primary charged particles at mid-rapidity is 
dh/drj = 1584 ± 4 ± 76 syst, an increase of about a fac- 
tor 2.2 to central Au-Au RHIC collisions at y^snn = 
0.2 TeV [1]. The study of particle yields per unit of ra- 
pidity obtained after integration of transverse momentum 
particle spectra eliminates the need to model the distor- 
tion of spectra introduced by explosive dynamics (see, 
e.g., [3]) of highly compressed matter created in high en- 
ergetic collisions. 

Following the reference [3], where deformation of the 
spectra were studied in model cases, we interpret the 
measured dh/drj to be equivalent to a central rapid- 
ity density dh/dy ~ 1800. This experimental result 
enables us to offer prediction for a variety of different 
particle multiplicities, dependent on adopted hadroniza- 
tion model and a first interpretation in terms of particle 
source bulk properties. 

SHM — Statistical hadronization model: QGP 
hadronic particle production yields are generally consid- 
ered within the statistical hadronization model (SHM) [3, 
1^ . SHM has been successful in describing hadron produc- 
tion in heavy ion collisions for different colliding systems 
and energies. Some view the SHM as a qualitative model 
and as such one is tempted to seek simplicity in an effort 
to obtain an estimate of the yields for all hadrons with 
just a small number of parameters 

Improving experimental precision, along with physics 
motivation based in qualitative dynamics of the 



hadronization process has stimulated refinements involv- 
ing a greater parameter set allowing to control the dy- 
namically established yield of different quark fiavors, gen- 
erally referred to as chemical non -equilibrium SHM 
This is achieved by introducing statistical occupancy pa- 
rameters 7i > 1, i — q, s, c, where s is the strange and c 
the charm quark flavor. It can be assumed that up and 
down quark yields q — u, d are equally equilibrated. We 
will not discuss the charm flavor here. 

Moreover, we are interested in precise description 
of the bulk properties of the particle source, such as 
size, energy and entropy content of the QGP fireball. 
This requires precise capability to extrapolate observed 
hadron yields to unobserved kinematic domains and un- 
observed particle types. This is the case for chemical non- 
equilibrium approach as demonstrated by the smooth 
systematic behavior of physical observables as a func- 
tion of collision conditions such as reaction energy Q or 
collision centrality 0. 

With increasing collision energy, the baryon content 
at central rapidity decreases rapidly. It is expected that 
there remains a small excess of matter over antimat- 
ter at central rapidity at LHC In SHM, this is 
governed by chemical parameters Ag, As (equivalent to 
Ms I Ms I the baryon and strangeness chemical potentials 
of matter). Determination of values of /is, /is must 
await additional experimental data. We present, for ori- 
entation of how the baryon-antibaryon yields vary, re- 
sults testing different values of the expected potentials 
choosing a fixed input value A^ = 1.0055, which results 
for in jXB ^ 2.9 ± 0.3 MeV, jis ^ 0.7 ± 0.3 MeV, and 
jjLB ^ 2.0 ± 0.2 MeV, jis - 0.45 ± 0.05 MeV in case of 
equilibrium and non-equilibrium models, respectively. 

Note that since the number of strange and antistrange 
quarks in hadrons has to be equal, A^, and thus /xg is de- 
termined in order to satisfy this constraint. Another con- 
straint is also implemented, the total charge per baryon 
has to be Q/B ~ 0.4, since the stopping of electrically 
charged matter (protons) within a rapidity interval is the 
same as that of all baryonic matter (protons and neu- 
trons) . This is achieved by a suitable tiny up-down quark 
asymmetry, a feature implemented in the SHARE suit of 
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programs we are using Q ■ 

In the usual procedure of statistical hadronization 
modeling, the particle yields are used to find best sta- 
tistical parameters. The SHARE suit of programs was 
written in a more flexible way to allow also mixed a fit, 
that is a fit where a few particle yields can be combined 
with the given 'measured' statistical parameters to ob- 
tain best fit of other statistical parameters. When we 
were developing SHARE, this feature was created since a 
parameter such as temperature could be measured using 
spectral shape and thus should be not fitted again in the 
yield description but be used as an experimental input. 
To be general, this feature was extended to all statistical 
parameters of the SHARE program. 

This feature allows us to perform a fit of the mix of 
statistical parameters and particle yields. Our procedure 
has been outlined before [15| and is, here, applied for the 
first time including LHC experimental data. This allows 
us to predict within the chemical equilibrium and non- 
equilibrium SHM models the differing pattern of particle 
production such that dh/dy is fixed. 

Importance of understanding chemical (non)- 
equilibrium: When chemical non-equilibrium is derived 
from the particle yields, generating in the fit to data 
7i ^ 1, this suggests a dynamical picture of an explosively 
expanding and potentially equilibrated QGP, decaying 
rapidly into free streaming hadrons. Without a signifi- 
cant re-equilibration, the (nearly) equilibrated QGP can- 
not produce chemically equilibrated hadron yield. The 
high intrinsic QGP entropy content explains why equili- 
brated QGP turns into chemically overpopulated (over- 
saturated) HG phase space. — The fast breakup of QGP 
means that the emerging hadrons do not have opportu- 
nity to re-establish chemical equilibrium in the HG phase. 

The differentiation of chemical equilibrium and non- 
equilibrium SHM models will be one of the challenges we 
address in our present discussion facing the SHM model 
and interpretation of the hadron production. One could 
think that resolution of this matter requires a good fit of 
SHM parameters to the data. However, with the large 
errors on particle yields and lack of sensitivity to 7^ , this 
is not easy. 

The light quark phase space occupancy parameter 7^ 
can be only measurable by determining overall baryon to 
meson yield, and this cannot be done without prior mea- 
surement of hadronization temperature T . When par- 
ticle yield data is not available to measure both T and 
7q, one can fit only ^s/lq to the data, which is then re- 
ported as 7s accompanied by the tacit assumption 7g = 1. 
Since 7^ (or ^s/lq) controls the overall (relative) yield of 
strange quarks, one expects and finds in most environ- 
ments 7s ^ 1 (or 7s/7i3 ^ 1) and a value which increases 
with system size, and often with energy. 

We recognize considerable physics implication of un- 
derstanding the value of 7^, as this relates directly to 
the measurement of T, related to the phase transforma- 
tion condition Ttr of QGP to hadrons, studied within lat- 
tice QCD. In the chemical equilibrium context, T ~ Ttr. 



On the other hand, the chemical non-equilibrium SHM 
implies rapid expansion and supercooled transformation 
and hence T < Ttr, with estimated difference at 10-15 
MeV [13. 

For lower heavy ion reaction energies as compared to 
LHC, one can also determine the baryochemical poten- 
tial at hadronization. This then allows to compare lattice 
QCD transformation condition, Ttr, /^B,tr [nlJilL with 



the data fit in the chemical equilibrium SHM in and 
chemical non-equilibrium SHM ^5|. One finds that lat- 
tice results are much flatter compared to the equilibrium 
SHM, that is Ttr drops off much slower with /LtB,tr- On 
the other hand, the non-equilibrium SHM parallels the 
lattice data 15 MeV below the transformation boundary. 
This favors, on theoretical grounds, the chemical non- 
equilibrium approach. 

LHC predictions assuming chemical equilib- 
rium: Just one observable, the number of charged par- 
ticles dh/dy ~ 1800, fixes within statistical model a 
large number of particle yields. This is allowing to test 
SHM model. We explore chemical equilibrium and non- 
equilibrium conditions in turn. We will show that many 
of the particle yields vary little as function of hadroniza- 
tion condition for a prescribed dh/dy. We present an 
unusual number of digits, a precision which has nothing 
to do with experiment but is needed to facilitate repro- 
ducibility of our results. We further state the propa- 
gation error of the error A dh/dy = 100 which is mostly 
confined to the volume, but is in some cases also visible in 
T and also further below in the chemical non-equilibrium 
parameters 7^, 7s, 

Our procedure for the equilibrium SHM is as fol- 
lows: in the chemical equilibrium model, see table HI 
we set 7g = 1, 7s = 1 and as noted also on other 
grounds Xq — 1.0055. As follows from the constraint of 
strangeness balance s — s — 0. All our results maintain a 
fixed Q/B = 0.4. We then see that, at fixed hadroniza- 
tion temperature T chosen in table U to be 159, 169, 
179, 189 MeV, the yield of charged hadrons dh/dy is, 
correlated to the source volume dV/dy. Volume varies 
strongly with temperature, since particle yields scale, for 
T ^ m, as VT^. Actually, since T ^ m condition is not 
satisfied, the volume is changing more rapidly so that 
T'^dV/dy = Const., k ~ 7.2. 

Since charged particle number is fixed, we also expect 
that entropy of the bulk is fixed and that is true; up 
to a small variation due to variation with T in relative 
yield of heavy hadrons, the entropy content is dS/dy = 
14800 ± 400. As temperature increases, the proportion 
of heavy mass charged particles increases, and thus the 
pion yield and even the kaon yield slightly decrease with 
increasing T. Baryon yields is most sensitive to T: VL 
doubles in yield in the temperature interval considered. 
Our choice of A, fixes for each hadronization temperature 
the per rapidity net baryon yield also shown in table HI 
— We believe our choice is reasonable and has been made 
also so that we see that there is no need to distinguish 
particles from antiparticles. 
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TABLE I: Chemical equilibrium particle yields at y^sjvjv ~ 
2.76 TeV. Top section: input properties; middle section: prop- 
erties of the fireball associated with central rapidity; bottom 
section: expected particle yields, and some select ratios. 
* signals an input value, and ** a result directly following from 
input value (combined often with a constraint). All yields, 
but TT^vD) without week decay feed to particle yields. Error 
in dV/dy corresponds to error in dh/dy equal 100. 





T* [MeV] 


159 


169 


179 


189 


7j 


1 


1 


1 


1 


7? 


1 


1 


1 


1 


A* 


1.0055 


1.0055 


1.0055 


1.0055 


10^{Xs - 1)** 


2.06 


1.45 


0.89 


0.39 


{Q/Br 


0.4000 


0.4000 


0.4000 


0.4000 


(s - s)* 


0.0000 


0.0000 


0.0000 


0.0000 


(dh/dy)* 


1800 


1800 


1800 


1800 


dV/dy [im^\ 


5285±147 


3452±96 


228d±64 


1538±43 


dS/dy 


15155 


14940 


14690 


14420 


s/ o 


n no/i c; 
U.Uz4D 


U.Uzoo 


U.UzOo 


U.Uz ; U 


[MeV/rm J 


64. 1 


100 


153 


231 


E/TS 


0.859 


0.86 


0.87 


0.87 


ID / TP 


0.164 


0.158 


0.153 


0.150 


-C// V [Ljev /im J 


U.oyz 


U.Doz 
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u.yy / 


1.04 


^- ^+ 
TT , 7r 




ooU 


oZl 


Q1 Q 
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K~ 


141.3 


140.8 


139.0 


136.8 




142.1 


141.6 


140.1 


137.8 


P 


53.6 


63.1 


72.0 


79.8 


P 


51.9 


61.2 


69.7 


77.3 


A 


30.0 


36.3 


42.1 


47.3 


A 


29.2 


35.4 


41.1 


46.2 




4.45 


5.47 


6.41 


7.23 




4.36 


5.38 


6.31 


7.14 


n- 


0.772 


1.038 


1.314 


1.586 


{B~B)** 


4.81 


5.60 


6.29 


6.88 


p 


92.4 


96.6 


99.1 


100.3 





19.0 


20.5 


21.4 


21.9 


K0*(892) 


42.6 


45.3 


46.9 


47.6 


K0*(892)/K- 


0.301 


0.322 


0.337 


0.348 


fli/K0*(892) 


0.446 


0.452 


0.456 


0.460 




942 


933 


925 


916 


■n 


110 


111 


111 


110 


V' 


9.67 


10.4 


10.8 


11.1 




1251 


1251 


1249 


1245 



Yield of strangeness is slightly increasing, this in- 
crease is in heavy mass strange baryons, e.g.. A, and 
this depletes slightly the yield of kaons. Overall the 
specific strangeness per entropy yield grows very slowly 
from s/S ~ 0.0245 to 0.027. Several ratios, such as 
0/K°*(892) ~ 0.45, where several effects compensate are 
nearly constant. 

We also show the post-weak decay tt^^ yield which 
is relatively large and independent of hadronization T. 
The decay 7r° — >■ 77 generates a strong electromagnetic 
energy component. 

LHC prediction within non-equilibrium SHM: 
Within the non-equilibrium hadronization approach, we 
need to further anchor the two quark pair abundance pa- 
rameters 7(j and 7^ . In absence of experimental data, we 
introduce additional hadronization conditions, the rela- 
tive strangeness yield s/S and hadronization pressure P. 
We vary P, see tableHH just as we varied T the hadroniza- 



TABLE II: Chemical non-equilibrium particle yields, each col- 
umn for different hadronization pressure. See caption of table 
Ufor further details. 



P* [MeV/fm^] 


60.3 


70.0 
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O.Ojoy 


0.0370 
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u.uuuu 
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u.uuuu 


nono 
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(dh /dv)* 

yu,iij u,y j 


1800 


1800 


1800 


1800 


1 [iViC V J 


±01 .z 


Lo^.o m u. 1 


10 1 . ( m u. ± 


loy .D m u. 1 


7? 


i .oyy 

±0.001 


1 finn 

i .DUU 

±0.008 


±.uu± 
±0.009 


i .oyy 
±0.011 


7s 


2.913 
±0.008 


2.842 
±0.030 


2.745 
±0.030 


2.721 
±0.016 


dV/dy [fm3] 


5469±542 


4731±136 


4043±119 


3705±168 


dS/dy 


13924 


13879 


13794 


13797 


E/TS 


1.060 


1.060 


1.059 


1.059 


P/E 


0.170 


0.168 


0.165 


0.164 


E/V [GeV/fm^] 


0.354 


0.417 


0.497 


0.550 


n~ , 7r+ 


858 


854 


850 


848 


K- 


192.0 


190.2 


186.5 


186.2 


K+ 


192.9 


191.2 


187.5 


187.3 


P 


32.9 


36.3 


40.3 


42.5 


P 


31.8 


35.2 


39.0 


41.2 


A 


28.9 


31.8 


34.8 


36.9 


A 


28.1 


31.0 


33.9 


35.9 




6.92 


7.56 


8.12 


8.60 




6.77 


7.40 


7.96 


8.43 


n- 


1.56 


1.73 


1.89 


2.03 


(B-B)** 


3.640 


3.973 


4.328 


4.539 


P 


56.1 


58.8 


61.7 


63.2 


4> 


30.0 


30.7 


30.8 


31.4 


K0*(892) 


39.9 


41.4 


42.6 


43.6 


K0*(892)/K- 


0.208 


0.218 


0.228 


0.234 


0/^0* (892) 


0.751 


0.741 


0.722 


0.721 


ttO 


988 


983 


979 


977 




134 


132 


128 


128 


n' 


10.4 


10.7 


10.8 


11.0 




1398 


1396 


1389 


1391 



tion condition in the chemical equilibrium model. In the 
chemical non-equilibrium SHM, two conditions suffice to 
narrow considerably the values of three SHM parameters 
(7q, 7s and T), but only if we insist that a third condition 
E/TS > 1 is qualitatively satisfied. 

Strangeness yield is a natural hadronization condition 
of QGP. We consider the ratio of strangeness per entropy 
s/S* in which coefficients and other systematic depen- 
dencies cancel. Since the entropy contents is directly re- 
lated to the particle multiplicity dh/dy, in our case s/S 
implies s-yield and directly relates to strangeness pair 
yield. In QGP, ss-pairs are produced predominantly in 
thermal gluon processes and their yield can be obtained 
within the QCD perturbative approach. In a study 
which was refined to agree with the strangeness yield 
observed at RHIC, we predicted the value s/S ~ 0.037 
for LHC [1^. We use here this result, noting that 
higher s/S values are possible, depending on LHC formed 
QGP dynamics. The QGP expected dynamic strangeness 
yield is considerably higher than the chemical equilibrium 
yield, table HI 0.0270 < s^q/S < 0.0245. The greater 
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strangeness content in QGP is, in fact, the reason be- 
hind the interest in strangeness as signature of QGP. 

The second condition arises from the observation that 
once the statistical parameters were fitted across di- 
verse reaction conditions at RHIC, the one constant 
outcome was that the hadronization pressure P = 82 
MeV/fm^ [l7|- Choice of pressure as a natural QGP 
hadronization constraint is further rooted in the obser- 
vation that the vacuum confinement phenomenon can 
be described within the qualitative MIT-bag model of 
hadrons introducing vacuum pressure is -Bmit = 58 
MeV/fm^, while in a bag- motivated fit to hadron spectra 
which allows additional flexibility in parameters one finds 
Sfit = 112 MeV/fm^ Clearly, a range of values is 

possible theoretically, with the hadronization condition 
P = 82 MeV/fm'^ right in the middle of this domain. We 
will use this 'critical pressure' hadronization condition as 
our constraint, but also vary it such that 60 < P < 90 
MeV/fm^ so that we can be sure that our prediction is 
not critically dependent on the empirical value. The pres- 
sure seen in equilibrium model case, table HI has a range 
64 < P < 230 MeV/fm3. 

The third constraint is not imposed in its precise value, 
but we require that hadronization occurs under the con- 
straint that E/TS > 1. In comparison, for the equi- 
librium case, table H we have 0.86 < E/TS < 0.87 a 
relatively small variation. The importance of this quan- 
tity E/TS > 1 as a diagnostic tool for explosive QGP 
outflow and hadronization was discussed in (l4| . In 
fact, we find that a reasonable and stable hadronization 
arises in chemical non-equilibrium within a narrow inter- 
val 1.059 < E/TS < 1.060. 

In table ini the outcome of this procedure is presented. 
We state the actual values of parameters for which so- 
lution of all constraints was numerically obtained, thus 
in first column pressure is not 60 but 60.3 MeV/fm'^. 
With rising hadronization pressure, the hadronization 
temperature rises, but it remains well below the phase 
transformation temperature. As we have discussed, the 
low value of T in the chemical non-equilibrium SHM is 
consistent with the dynamics of the expansion, the flow 
of matter reduces the phase balance T. This in turn is 
then requiring that the light quark abundance parame- 
ter 7, ~ 1.6. This is the key distinction of the chemical 
non-equilibrium. It further signals enhancement of pro- 
duction of baryons over mesons by just this factor. Note 
that Js/lq — 1.72. This large ratio means that the yield 
of A and p do not differ much. This indicates strong en- 
hancement of strangeness, a first day observable of QGP 
formation, along with enhancement p^ . 

Comparison of SHM results: The large bulk 
hadronization volume dV/dy ~ 4500 fm"^ is suggesting 
that there will be noticeable changes in the HBT ob- 
servables allowing to produce such a great hadronization 
volume. The bulk energy content is found in both ap- 
proaches to be dE/dy = E/V x dV/dy = 2.00 ± 0.05 
TeV per unit of rapidity at y = 0. This is the thermal 
energy of QGP prior to hadronization measured in the 



local fluid element rest frame. 

The entropy content in the bulk for non-equilibrium, 
dS/dy — 13860 ± 64, is 5% smaller compared to equilib- 
rium case. This is due to the fact that non-equilibrium 
particle yields do not maximize entropy. We note that the 
yields of pions, kaons, and even single strange hyperons 
are remarkably independent of hadronization pressure, 
or, in the equilibrium case, temperature, even though the 
volume parameter changes greatly. This effect is coun- 
teracted by a balancing change in hadronization temper- 
ature since the yield of charged hadrons is fixed. The 
hadronization energy density is very close to E/V ~ 0.5 
GeV/fm'^, and it tracks the pressure since the ratio P/E 
is found to be rather constant. 

The yield of multistrange hadrons is much enhanced in 
chemical non-equilibrium model, compared to the equi- 
librium model, on account of 50% increased yield of 
strangeness, which is potentiated for multistrange par- 
ticles as was predicted to be the signature of QGP 19} . 
The yield of A, for the most favored hadronization condi- 
tion in both equilibrium and non-equilibrium, can in fact 
be lower in the non-equilibrium case than in the equilib- 
rium, yet at LHC the yield of K is always 40% greater. 
Like in the equilibrium results, we observe that several 
yields are largely independent of the hadronization condi- 
tion, meaning that ratios such as (f>/h could be a distinc- 
tive signature of hadronization, differentiating the two 
primary models. This is illustrated in figure [T] where the 
ratio is shown as function of resultant hadronization T. 
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FIG. 1: (color on-line) Specific yield </i//i for chemical equi- 
librium model to right bottom (blue) and chemical non- 
equifibrium model left up (red). 

One can, however, argue that (p/h, seen in Fig.[Tl could 
be brought about by strangeness enhancement, not re- 
quiring that 7^ > 1, for a more complete discussion 
see [20|. To narrow the choices, we propose to study 
two more ratios, shown in Fig.[2j The top section shows 
K*/K~, as function of the yield K~ /h. The lower right 
(red) non-equilibrium result shows strangeness enhance- 
ment at low hadronization T since K*/K~ mainly de- 
pends on T. We have shown both particle and antipar- 
ticle ratios derived from our fixed input for net baryon 
yield to illustrate that differentiation of these results will 
not be easily possible. 
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FIG. 2: (color on-line) Top: ratio of resonance K* to kaon 
yield K, as function of specific kaon yield K~ /h, the equi- 
librium model (blue) is left up and chemical non-equilibrium 
(red) is bottom right. Bottom frame: the K/tt ratio as func- 
tion of p/K ratio. Bottom right (blue) is equilibrium model 
and upper left is non-equilibrium model. 



The bottom section, in Fig. [5J shows the strangeness 
enhancement in format of K/tt as function of the easiest 
to measure baryon to meson ratio which is a measure of 



absolute magnitude of 7, here specifically, 7g/7s7g- We 
see the equilibrium model to lower right (blue) while the 
non-equilibrium model is upper left. 

Summary and conclusions: We have obtained par- 
ticle yields within the statistical hadronization model for 
the LHC-ion run at ^/snn = 2.76 TeV. We have discussed 
both bulk properties of QGP at breakup and the resulting 
particle yields. These can vary significantly depending on 
the hadronization mechanism. Distinctive features asso- 
ciated with QGP based strangeness enhancement and fi- 
nal state chemical non-equilibrium were described and 
strategies leading to better understanding of chemical 
(non-)equilibrium were proposed. 

We have shown how enhanced yield of strangeness with 
the phase space occupancy 7^ ~ 2.75 when 7g ~ 1.6 mod- 
ifies the yield of strange hadrons and detailed predictions 
for the observables such as cj)/h, K*/K, K/tt, p/K, A/p 
were offered. Enhanced yields of (multi)strange particles 
are tabulated. We note that absolute yield of cj) is en- 
hanced by a factor 1.5 in the non-equilibrium compared 
to equilibrium hadronization. There is no significant de- 
pendence of the (j) yield on hadronization condition mak- 
ing it an ideal first day differentiating chemical equilib- 
rium from non-equilibrium. 

The large bulk hadronization volume dV/dy ~ 
4500 fm'^ related to HBT observables, the local rest frame 
thermal energy content dE/dy\o = 2 TeV constrains hy- 
drodynamic models. A large yield of tt'^, rj and thus of 
associated decay photons is noted, enhanced somewhat 
in the. chemical non-equilibrium case. 
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